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SUMMARY 

/ 

The governing equations for an arbitrary rigid body sliding on a landing 
surface are used to derive a stability criterion which relates the critical, 
values of initial velocities to the coefficient of friction, center-of-gravity 
location, and initial angle of contact. A numerical application of the sta- 
bility criterion is made for a vehicle used in an experimental investigation. 


INTRODUCTION 


A critical problem area in the design of a space vehicle is the landing 
system. As conventional landing gears would be quite heavy, new methods for 
absorbing landing impacts have to be devised. One technique that contains 
little weight chargeable to the landing system is the one referred to as the 
skid-rocker landing system. (See, for example, refs. 1 and 2.) In this sys- 
tem the vehicle is permitted to slide on its curved lower surface, and the 
friction between this surface and the landing surface is relied on to stop the 
vehicle eventually. 

An application of this landing concept is discussed in reference 1. The 
skid-rocker landing characteristics of two proposed space vehicles with dif- 
ferent un dersurfaces were evaluated experimentally through the use of dynamic- 
model testing techniques. The results of this investigation showed that upon 
initial contact with the landing surface, the friction force at the contact 
point could couple with the inertia forces in such a way as to cause the vehi- 
cle to tumble. As such an instability could preclude the use of a specific 
configuration, the stability characteristics of a landing vehicle become an 
important design consideration. 

It would appear to be unrealistic to consider determining experimentally 
the stability characteristics of all competing vehicle shapes since the number 
of such vehicles is great and each vehicle should be optimized with. respect to 
variations in the shape of the sliding surface and its weight distribution. 
Hence some analytical method of determining the stability of sliding bodies is 
needed. The presentation of such an analysis is the purpose of this paper. 


This paper is divided into two parts. In the first part the governing 
equations for an arbitrary rigid body sliding on a surface are presented, along 
h the initial conditions both before and after contact. These nonlinear dif- 
erential equations are solved exactly to obtain closed-form expressions for the 
angular position in terms of the initial conditions. These expressions in turn 
are used to derive the stability criterion. The second portion of the paper 
consis s of numerical results obtained by the application of the derived sta- 
rt^ 7 pterion to on e of the vehicles used in the experimental investigation 
iscussed in reference 1. Numerical results showing the influence of coeffi- 

of e the 0 vph^'r'^ 10n, ^ nklng ® pee ^ and cen- ter-of -gravity location on the stability 
vehicle are discussed. In addition, comparisons between theory and 
experiment are presented. J 


SYMBOLS 


b 


perpendicular distance from body reference line to 
(fig. 1) ” 


center of gravity 


b = b/D 

U distance from base of bo <iy reference line to center of gravity 

measured parallel to body reference line (fig. l) 

c" = c/D 

D base diameter of vehicle (fig. l) 

d horizontal distance from contact point to center of gravity (fig. l) 

F friction force 

g acceleration due to gravity 

h vertical distance from contact point to center of gravity (fig. l) 

value of h at initial contact 
I polar mass moment of inertia 

N normal force 

P,Q functions of 9 defined in equations (10) and (ll) 

t time 

Vh horizontal velocity 

Vv vertical velocity (positive downward) 
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v e 

w 

x 

x o 

Z>l 

a,T) 

9 


0L 


rotational velocity 
weight of vehicle 

coordinate of center of gravity, measured parallel to landing surface 
value of* x at initial contact 
coordinates of undersurface 
dummy variables 

attitude of vehicle (angle between body reference line and the 
perpendicular at the point of contact, fig« l) 

angular limit of stability 


0 O initial-contact angle 

^ friction coefficient 

fl = wr 

co angular velocity, 9 

cdq value of tn at initial contact 

co t critical value of to at initial contact 

O, 

Dots on symbols denote differentiation with respect to t; primes denote 
differentiation with respect to 0 . 


ANALYSIS 


Description of Problem 

The problem to be considered is that of the motion of a rigid body sliding 
on a landing surface in such a way that continuous contact is maintained between 
the body and the landing surface. Only bodies having sliding surfaces with no 
reverse curvature are considered; thus only one point is in contact with the 
landing surface at any time. The body is assumed to have no roll or yaw upon 
contact, so that all forces including the inertia forces are coplanar . 

A typical cross section through the center of gravity of such a body is 
shown in figure 1. The position of the body at any time t can be given by the 
coordinates x and 0, where x is the coordinate of the center of gravity . 
parallel to the landing surface and 0 is the angle between some reference line 
in the body and the perpendicular at the point of contact. The position of the 
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center of gravity with reference to the point of contact is given by h, the ver- 
tical distance above the surface (positive up)^ and d, the horizontal distance 
from the contact point (positive aft) . Note that for a given shape of undersur- 
face both h and d are known functions of 0 only. The location of the center 
of gravity with reference to the base of the body reference line is given by the 
distances b and c, measured perpendicular and parallel to the body reference 
line respectively. The quantities b, c , d, h, and 0 are shown in the posi- 
tive sense in figure 1. Note that the contact angle 0 is the negative of the 
conventional angle of attack in aerodynamics. 



The forces acting on the body consist of the normal force N, the friction 
or drag force F, and the weight of the body W. All forces are positive in the 
directions shown in figure 1. 



Governing Equations 


Application of Newton's laws of motion results in the following equations: 


= -Sfe 

(1) 

g 


Nd = 19 

( 2 ) 

W = *h 

( 3 ) 


g 


where g is the acceleration due to gravity and I is the mass moment of 
inertia about the axis through the center of gravity and perpendicular to the 
plane of motion. The dots over the symbols refer to differentiation with respect 
to time t. 

The drag force F is related to the normal force by the equation 


F = pN 


(*0 


where p is the coefficient of friction, which is assumed to be constant. Note 
that since h is a function of 8 alone, and 9 is a function of t, the 
vertical velocity and acceleration of the center of gravity are given by 


h = h'0 ( 5 ) 

and 

h = h"0 2 + h’0 (6) 

where the primes denote differentiation with respect to 9 . 

Angular Motion of Vehicle 

Solution of equation (2) after substitution of equation ( 4 ) results in an 
expression for the normal force in terms of 9: 

N = - 0 (T) 

d - ph 

With this expression for N and equation ( 6 ), equation ( 3 ) can written as 
follows : 
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-w 


( 8 ) 



+ + ^"e 2 = 

g / g 


Since the quantities d and h are functions of 9 and the q uan tities I, 
W, hj and g are constants, equation (8) constitutes the sole governing equa- 
tion for 9. The solution of this equation gives the orientation of the body at 
any time. Note that equation (l), which involves the acceleration in the hori- 
zontal direction, does not enter into the derivation of equation (8) • Hence the 
angular motions of the vehicle are entirely uncoupled from the horizontal motion. 
The horizontal velocity, however, must at all times be high enough to insure that 
the resultant velocity of the contact point is positive so that the friction 
force does not change direction. 

Before proceeding with the solution of equation (8), it is convenient to 
rewrite the equation as 


29 + P9 2 = Q 


(9) 


where both P and Q are functions of 0 defined by 


P = 


2 | 1 " 


1 + *h’ 


d - |ih g 


( 10 ) 


Q = - 


2W 


-I— + 


d - jih S 


(id 


Although 0, the angular coordinate, is necessary in obtaining the orienta- 
tion of the body, a more important variable for the stability problem is the 
angular velocity u> where 


9 = 0 ) 


and 


0 = (h = OKI} 1 

With equations (l2) and (13), equation (9) can be written in terms of u>: 

20*0)' + p = Q 


(12) 

(13) 

OM 


This nonlinear equation (eq. (l4)) can be reduced to a linear differential equa- 
tion by the transformation 
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( 15 ) 


Q. ^ af 


so that 


fl’ = 2<aoi' 


( 16 ) 


Substitution of equations ( 15 ) and (l6) into equation (l^) results in the 
following linear differential equation in fl: 


fl' + Pfl = Q 


( 17 ) 


With the use of the integrating factor 
tion ( 17 ) can be written as 



the solution to equa- 



fle 


- fl( o) 




or, in terms of to, 


2 u 
(xre 


P(a) da 


0) 


-( 0 ) 



Q(a)e 


/ 0 ap (i )d i 

s da 


(18) 


The value of 
conditions at t = 


0 )^( 0 ) is obtained by evaluating equation (l8) for the initial 
0; that is, for 9 = 0 Q and o>(9 0 ) = <u 0 : 


a£(o) 



o 

P(a)da 



Q(a)e 


fo 

> da 


( 19 ) 


With the use of equation ( 19 ), equation (l8) can be solved for o^: 



Equation (20) can be used to obtain the value of the angular velocity <x> 
for any value of 9 in terms of known initial values of 0 O and o> 0 . The 
angular acceleration can be obtained by differentiating equation (20) with 
respect to 9 and making use of equation (13) • Thus 
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0 =-|p(9)e 


P( Tl)d-Tl 


f P(T))dT] p 9 
03 Q 2 e + / Q(a)e" 

JQr, 


f n PUUn 

00 A 

da 


+ |Q(0) (21) 


or 


9 = - — Rjd 2 + 

2 2^ 


( 22 ) 


Equation (22) could also have been obtained directly from equation ( 9 ) by noting 
the relationship in equation (l2) . 


Stability Criterion 

Upon initial contact of the vehicle with the landing surface, friction force 
and normal force are generated at the contact point. These forces couple with 
the vertical inertia forces to create angular accelerations of the vehicle. If 
the initial conditions and geometry of the vehicle are such that these accelera- 
tions are large enough, the vehicle will topple over. If these forces are not 
large enough, the vehicle will reach a maximum angle of tip and then oscillate 
about some mean position until the friction forces bring it to a halt. 

From this description of the vehicle motion it can be seen that for a 
vehicle to be stable two conditions must be met. First, the vehicle must oscil- 
late between two extreme values of 8; consequently, there must be two values of 
9 at which the angular velocity vanishes: 

9 (^E+) = m (^E+) = 0 ( 23 a) 

9 ( 0 E-) = fD ( 9 E-) = 0 (25b) 

Second, the angular acceleration at these points must be such as to cause the 
vehicle to return to the mean position; that is, 

9 ( e E+) < 0 (24a) 

e(e E _) > 0 ( 24 b) 


In these expressions 0^+ is the value of the extreme on the positive side of 
the mean position and 0 E _ the value of the extreme on the negative side. The 
angular limit of stability 9^ is obtained when at one of these extremes the 
angular acceleration attains a value of zero. Thus 0^ is given by 
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e(e L ) = ai(e L ) = o 


( 25 a) 



( 25 b) 


For each vehicle considered, an investigation must he made of the accelerations 
at the extreme angular positions in order to determine which extreme yields the 
limit of stability. 


Equation (25a) , upon substitution from equation (20), yields 




P ( T] ) dTj 


?( T]) d_Tj 


la 


(A 


>o,L 


+ e 


Q(a)( 


1 0 


P(n)dT 


da 


= 0 


from which the critical initial angular velocity fD o,L can determined. 

fe l 

p(n)dn 

t o 


Since 


U Q 

e °o 


the equation for becomes 


%,L = - e 


) oa 

P(r 1 )d T i TSl L PUAn 

/ Q(a)e^° da 

JSn 


(26) 


In addition, as 0 ^ must satisfy equations (25a) and (25b) simultaneously 
(see eq. ( 22 )), 


Q(e L ) = 0 (27) 

From the definition of Q in equation (ll) it can he seen that equation (27) is 
satisfied only if either 


h’ = oo 


or 



(28) 


The former condition is impossible when it is noted that 

h' = d 


(29) 


9 



Consequently the latter (eq. (28)) must be one of the governing stability 
requirements . 

Note that when d and h are related as in equation (28), the resulting 
force at the point of contact intersects the center of gravity of the body. For 
most practical shapes of undersurfaces used as skid rockers there are three solu- 
tions to equation (28) (that is, three positions at which eq. (28) is satisfied). 
The intermediate solution is associated with a position that is always stable 
while the extreme values are the possible limits of stability. 

The attitude of a given vehicle at the limit of stability is obtained from 
equation (28) and is a function of the coefficient of friction and the shape of 
the body. Once this limiting angle is known, the critical initial angular 
velocity to 0 ^ is found as a function of initial contact angle 0 Q from equa- 
tion (26) . The integrands appearing in this equation, however, are relatively 
complex functions of the shape in terms of h, h', and h" and consequently 
must be integrated numerically. 


Initial Conditions 

The initial conditions to be used in equation (26) are, of course, the con- 
ditions just after contact is made with the landing surface. In many instances, 
however, the values of the velocities which are known are those just prior to 
contact. These values are related to the initial conditions (just after contact) 
through impulse momentum equations: 





N dt 




(50) 



- a)at = i(a> 0 - v 0 ) 


(3i) 



w)at = |(h 0 + v v ) 


(32) 


where V^, V y , and Vq are, respectively, the horizontal, vertical, and rota- 
tional velocities just prior to contact, and x Q , h Q , and u> 0 are the corre- 
sponding initial velocities (just after contact). The integrals on the left of 
equations (30) to (32) represent the impulse existing between the landing sur- 
face and the vehicle. Their actual value depends on the resilient properties of 
both the vehicle and the landing surfaces. However, the time interval At, 
representing the duration of the impulse, is assumed to be small enough so that 
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rat r At 

/ N(qh - d)dt = (qh - d) / N dt 

0 Jo 


(53a) 


rat Pat 

W dt = W dt = 0 

'0 Jo 


(33b) 


The initial vertical and rotational velocities are related through equa- 
tion ( 5 ): 

bo = b ' ^9 G ) «*> ( 3^ ) 

With the use of equations (33b) and ( 3 ^-), equation (32) reduces to 


’'at 


N dt = — 


< 0 


( 9o) fJi o V v 


(35) 


The initial angular velocity tt) 0 , in terms of the velocities just prior to 

contact, is obtained from equation ( 31 ) after making use of equations' (33a) and 

(35): 


(JLU 


-_d( 

e 0 ) 

1 - q h| 

( 9 °) 

J J 

V + Wq 

I + 

CD 

1 

^5^ 

Bo) 

OQJ^ 

CD 

O 


(36) 


A similar expression for x 0 is obtained with the use of equations ( 30 ), (35 ) > 
and ( 36 ) : 


x o = v h “ ^ 



Vy 

+ h'(e 0 

he 

I + 

d| 


1 - H h( 

0 o) 

*h'l 

g 

w 


(37) 


Although the values of the impulse quantities depend on the resilient prop- 
erties of the contacting surfaces, the initial velocities just after contact can 
be obtained in terms of the velocities just prior to contact without knowing 
these resilient properties. 


Solution Neglecting Friction 

As was previously stated, the integrals (see eq. ( 26 )) that must be evalu- 
ated to obtain the critical initial angular velocity are integrable only for 
vehicles with specific simplified sliding surfaces. If the effects of friction 
are neglected, however, the quantities inside the integral signs become exact 
differentials; consequently, the integrals can be evaluated exactly for the 
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general case. As the stability criterion is obtained in a simple closed-form 
expression, the solution will be presented herein. 

As seen from equation (28), the attitude at which instability occurs when 
friction is neglected, 0p, is the position at which the center of gravity is 
directly above the point of contact - that is, the attitude for which d, and 
consequently h 1 , is zero. With the attitude of the vehicle at the point of 
instability known, the critical initial angular velocity at which the instability 
occurs, I s obtained from equation (26). 

With the assumption that the friction coefficient \i is zero, the expres- 
sions for P and Q (eqs. (10) and (ll)) reduce to 


2 Wh’h” 

Ig + w(h ’) 2 


and 


Q = -gWgfr 1 

Ig + w(h ') 2 

Note that P is now an exact derivative; therefore 


(38) 


(39) 



P d-q = log e U[ig + W(h')d 

^O 


( 40 ) 


where 

R c = Ig + w[h'(o)] 2 

Also 



f[ig + w(h') 2 ] 


OD 


(te) 


By the use of equations ( 39 ) and ( 4 - 2 ), equation (26) can be integrated to 
yield 


2 = 2Wg[h(e L ) - h(e Q j] 
°’ L ig + w[h’(e 0 )] 2 


(^ 3 ) 


12 


Consequently, all values of initial velocities before contact that result in 
initial angular velocities cjd 0 (as calculated by eq. (56) with p = 0) such that 
(jo 0 ^ > 03 o result in instability. 


NUMERICAL RESULTS 

The purpose of this section is twofold: first, to illustrate the applica- 

tion of the stability criterion In obtaining a stability boundary for a typical 
sliding body, and second, to show the quantitative effects of various parameters 
on the stability limit of a typical vehicle. 


Description of Vehicle 

The cross section of the vehicle whose stability is to be investigated in 
this section is shown in figure 2. The shape was taken to be the same as that of 



Figure 2.- Cross section of configuration C of reference 1. 

configuration C used in the experimental investigation discussed in reference 1. 
The undersurface is symmetrical and consists of hyperbolic and circular segments 
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having a common tangent at their junctions. The equations for the curves and the 
expressions for h in the three segments of the undersurface shown in figure 2 
are as follows: 


For the hyperbola (-22— < 0 < 22^ ): 


= -1.784 + J 0 . 23051 2 + 1.784 2 


h = (c + 1.784) cos 9 - b sin 0 - 1.784^cos 2 9 - 4.3583 sin 2 0 


~i o 

For the circle at 0 > 22^ : 

3 


Z = 1.977 - {( 0.3533) 2 - (t - 6.052) 2 


h = (c - I. 977 ) cos 0 - (b - 6 . 032 ) sin 9 + 0.3353 

I 0 

For the circle at 0 < - 22— : 

3 

Z = 1.977 - \j(o. 3333) 2 - (0 + 6.032) 2 

h = (c - I. 977 ) cos 9 - (b + 6 . 032 ) sin 0 + 0.3335 

The value of the ratio of the mass to the mass moment of inertia is 

— = 0 . 09 ^ per sq ft 
Ig 


Application of Stability Criterion 

The first step in obtaining a stability boundary is to calculate the 
angular limit of stability 0 ^ from equation ( 28 ) for the assumed value of the 
coefficient of friction. As the vehicle in figure 2 is symmetrical, it is con- 
venient to use the line of symmetry as the angular reference line. For this 
symmetrical body, with the center-of-gravity positions considered, the highest 
of the three values of 0 ^ from equation ( 28 ) yields the limit of stability. 
Values of 0 l for various values of the coefficient of friction are shown in 
figure 3 . With the value of 0^ known, (the critical initial angular 

velocity after contact for which instability or tumbling occurs) is calculated 
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c - . 24 

*B-o 


for a given contact angle 0 O from equa- 
tion ( 26 ). The integrals appearing in 
equation (26) are evaluated numerically 
by use of trapezoidal integration with 
2 / 3° intervals of 0. 

Any combination of sinking velocity 
V y and of pitching velocity Vq that 

results in a value of ad 0 ^ (as calculated 
from eq. (56)) greater than u> Q ^ leads 

to instability. As the experimental tech- 
nique used in reference 1 does not permit 
the vehicle to have a pitching velocity, 

Vq is assumed to be zero in the numerical 

analysis used herein. Consequently, the 
value of u> 0 in the following discussion 

will be a function only of the sinking 
velocity V v and the attitude at impact 

as given by 0 Q . 



I 1 1 1 i 1 1 i 1 

0 .2 .4 .6 .8 

H- 

Figure 3 *- Variation of angular limit of 
stability Gl with friction coeffi- 
cient 1 1 . 


Discussion of Results 

Stability boundaries were calculated for the vehicle of figure 2 for a 
range of values of sinking velocity, coefficient of friction, center-of -gravity 
location, and contact attitude. The results of these calculations are shown in 
figures 4 to 7 - 

In figure 4 the critical stability curve is shown as a function of sinking 
velocity and contact angle for a vehicle with given center-of-gravity location 
and coefficient of friction. For all combinations of contact angle and velocity 
that fall below this curve the vehicle is stable. For values that fall above 
this curve the vehicle tumbles. The distance between the two curves for a given 
velocity represents the range of contact angles for which a stable landing is 
possible, and is referred to as the region of safe touchdown. 

The amount of attitude control necessary for the landing vehicle decreases 
with an increase in the size of this region of safe touchdown. For example, at 
a sinking velocity of 40 feet per second, the controls must be sufficiently 
accurate to land the vehicle within a range of 10 ° of contact angle. At a veloc- 
ity of 20 feet per second this range increases to 20° and thus less control is 
necessary. For velocities below 20 feet per second there is a sudden increase in 
the size of the safe touchdown region. 

Note that for sinking velocities greater than 20 feet per second some kind 
of attitude control is always necessary. As can be seen in figure 4 , for veloc- 
ities above 20 feet per second the only stable contact angles are positive; con- 
sequently the landing of the vehicle must be controlled to assure such an angle. 
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Even at lower sinking velocities the 
desirability of some kind of attitude 
control is indicated by the fact that the 
range of positive safe contact angles is 
considerably larger than the range of neg- 
ative angles. 

The influence of the coefficient of 
friction and the height of the center-of- 
gravity location (c in fig. 2) on the 
region of safe touchdown is shown in fig- 
ures 5 and- 6, respectively. In figure 5 
the critical stability curve of the 
vehicle is shown as a function of the 
contact angle and coefficient of friction 
for a sinking velocity of 10 feet per sec 
ond and the same center-of-gravity loca- 
tion used to obtain the results of fig- 
ure 4. In figure 6 the stability curve 
is shown as a function of the contact 
angle and a dimensionless center-of- 
gravity height c = c/D, for the sinking 
velocity of 10 feet per second (same as 
used in fig. 5) and coefficient of fric- 
tion of 0.40 (same as in fig. 4). 




Figure 4.- Influence of sinking 
velocity V v on stability. 



Figure 5*- Influence of friction 
coefficient ^ on stability. 


Figure 6.- Influence of center-of-gravity 
height 'c on stability. 
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From figures 5 and 6 it can be seen that, the region of safe touchdown 
becomes larger with a decrease in both the coefficient of friction and the height 
of the center of gravity. Furthermore, in both figures, for reasonable values of 
li and cT the range of positive safe contact angles is considerably larger than 
the range of negative angles. This fact supports the desirability of attitude 
control. 


The influence of offsetting the center of gravity from the line of symmetry 
is shown in figure 7. In this figure the critical stability curve is given as a 
function of the contact angle and a dimensionless offset distance b = b/D. The 
values of c - and jjl are the same as in figure 4 and V v is equal to 10 feet 

per second. This figure indicates that, for a sufficiently high negative offset, 
the range of negative safe contact angles might be made large enough to eliminate 
the necessity for attitude control. This situation, however, would require that 
the location of the center of gravity with respect to the motion of the vehicle 
be made a design criterion. 


For completeness the comparison between theory and experiment presented in 
reference 1 is shown in figure 8. In this figure the curves show the calculated 
stability boundaries as functions of center-of-gravity height and landing atti- 
tude for three values of friction coefficient and a sinking velocity of 10 feet 
per second. The symbols show the results of the experimental investigation; the 
open symbols represent stable landing attitudes while the closed symbols repre- 
sent unstable ones. In general the agreement between theory and experiment is 
good. 



% 

Figure 7.- Influence of center-of-gravity 
offset b on stability. 



Figure 8.- Comparison of theoretical results 
with experimental results for configura- 
tion C of reference 1 showing the influence 
of center-of-gravity height ~c on stabil- 
ity. Vertical velocity V v = 10 ft/sec; 

center-of-gravity offset b = 0. 
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CONCLUDING REMARKS 


A method of calculating a stability criterion for sliding bodies has been 
developed. The method resulted from an exact solution of the nonlinear equations 
of an arbitrary rigid body sliding on a surface in such a way that continuous 
point contact is maintained. The range of values of initial contact angles which 
result in stable motion is obtained from the solution of two equations, one of 
vhich involves integrals that, for most cases, must be integrated numerically. 

The method was used to obtain the stability boundaries of a vehicle used in 
an experimental investigation. From the numerical results it was concluded that 
although the region of safe touchdown could be increased by proper choice of 
center-of- gravity location and sinking velocity, some attitude control for the 
landing vehicle is desirable. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va. , January 15, 1963 . 
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